Storage of spin squeezing in a two-component Bose-Einstein condensate 
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Efficient control of spin squeezing in a two-component Bose-Einstein Condensate is studied by 
rapidly turning-ofT the external field at a time that maximal spin squeezing appears. We show 
that strong reduction of spin fluctuation can be maintained in a nearly fixed direction for a long 
time. We explain the underlying physics unambiguously, and present analytical expressions of the 
maximal-squeezing time. 
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Spin squeezing has attracted much attention for 
decades not only because of fundamental physical in- 
terests 0, H, H, 0, Hi, but also for its possible appli- 
cation in atomic clocks for reducing quantum noise @ 
and quantum information @, 0, H, S G3 ■ Formally, 
the spin squeezing is quantified via a parameter £ = 
(AJ n ) m in/-\/j/2, where (AJ„) m i„ represents the small- 
est variance of a spin component J n = J • n normal to 
the mean spin (J). For the coherent spin state (CSS), 
the variance (AJ n ) m i n = y/j/2 (i.e., £ = 1). A state is 
called spin squeezed state (SSS) if its variance of the spin 
component is smaller than that of the CSS, i.e. £ < 1. 

Kitagawa and Uea have investigated the spin squeezing 
generated by the so-called one-axis twisting (OAT) model 
with Hamiltonian: -ffoAT = 2k J 2 Possible realiza- 
tion of the OAT-type spin squeezing in a two-component 
Bose-Einstein Condensate (TBEC) has been proposed in 
Refs. [1, S0rensen et al. also argued that macro- 
scopic quantum entanglement can be characterized by 
using the OAT-type spin squeezing in the TBEC i. Be- 
sides the TBEC, most recently, Takeuchi et al. consid- 
ered another realization of the OAT-type spin squeezing 
by using the interactions between atoms and off-resonant 
light (paramagnetic Faraday rotation) [l^]. To coher- 
ently control spin squeezing, Law et al. introduced ad- 
ditional Josephson-like (or Raman) coupling to the OAT 
model: Hr — 2k J 2 + £IrJ x [l3(- Such a model has also 
been used to prepare arbitrary Dicke states [l4T |. 

In addition to the generation of the SSS itself, it is de- 
sirable to maintain not only the squeezing but also its di- 
rection for a long time Q ■ Jaksch et al. have shown that 
the OAT-type SSS can be stored for arbitrarily long time 
by removing the self- interaction [la ]. However, it might 
not be easy to handle in experiment since the precisely 
designed additional pulses are crucially required. In this 
letter, we propose a simple mechanism to obtain long- 
lasting spin squeezing in the TBEC. Our scheme is quite 
easy to realize in experiment since it can be achieved 
by turning-off the Josephson coupling once the TBEC 
reaches its maximal spin squeezing. 

We consider a two-component weakly-interacting BEC 
[l6l . Kl\ consisting of N atoms in different atomic hy- 
perfine states |o) and \b) coupled by a time- varying mi- 



crowave field with Rabi freq uency fl r f- Based on the 
two-mode approximation jl8l . [T9L I20I l2ll I22I |23| , the to- 
tal Hamiltonian can be described by (h = 1): 



H(t) = 2kJ* + fl(t)J x , 



(1) 



initial CSS \j,-j) x = 
amplitudes c m (0) = — 



~3) 



then the initial 
Since the initial 



where k = {n aa + Kbb - 2k &)/4, and n a p = 
g a p J d 3 r \<j> a (r)<j)p(v)\ , with g a(3 = 47ra Q/3 /m (a, (3 = 
a, b) being the s-wave scattering strengthes between 
atoms. The condensate-mode functions 4> a normalized 
to unity satisfies a coupled Gross-Pitaevskii equations 
[24[. Here we focus on the case that the external field 
is turned off rapidly at a certain time to, so the time- 
dependent Josephson-like coupling can be written as 
Q(t) = tt R O(t - i), where Q R = Q rf J d 3 r<t)* a (r)(j) b (r) 
and O(i) is the step function. 

The state vector at arbitrary time t can be expanded 
in terms of eigenstates of J z : \ip(t)) — J2 m Cm(t) \ j,m), 
where —j < m < j and j = N/2. The equations of 
motion for the amplitudes c m (t) are obtained by solving 
time-dependent Schrodinger equation. We consider an 
= e- inJ y/ 2 \j, 

-iy+™ ( 2j 
23 \ j + m 

CSS satisfies c_ m (0) = c m (0) for even N, and c_ m (0) = 
— c m (0) for odd N, we can prove the mean spin always 
along the x direction. In addition, we will consider only 
positive k case by assuming a aa + abb > 2a a &. However, 
our results keep valid in the opposite case by using initial 
maximum weight state of J x , i.e., \j, j) x . 

Now let us briefly explain the basic principle of our 
scheme. The initial CSS can be prepared by applying a 
short 7r/2 pulse to a single-component BEC with all the 
atoms being in the internal state \a) [1, [23j]. After that, 
the external Josephson field is immediately switched on, 
so dynamical evolution of the spin system is governed by 
the Hamiltonian |1| with fl(t) — fl R . If the coupling is 
optimally chosen, the Josephson interaction results in an 
enhanced spin squeezing compared with that of the OAT 
For N = 10 3 , we find that the maximal squeezing 
£0 = 8.7076 x 1CP 2 can be obtained by choosing arbitrary 
f?ii in a region 10.777 < £Ir/k < 10.818. As shown by 
the dashed lines of Fig. [TJ the squeezing £ and the mean 
spin (J x ) show collapsed oscillations [25]. At the time 
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to, £ decreases to its local minimum £o with m j n = 0, 
while (J x ) increases to its maximal value (J x )o- The 
basic features of our scheme are exhibited by the solid 
lines of Fig. [T] We find that if we turn off the Jesephson 
field at the time to, the maximal squeezing £o can be 
stored in a fixed direction (i.e., (9 m m = 0) for a long time. 
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FIG. 1: Time evolution of (a) the squeezing parameter, (b) 
the squeezing angle, and (c) the mean spin (J x ) for N — 10 3 
and Q,r — 10.8k. Dashed lines: constant-coupling case; Solid 
lines: turning-off external field at the time Kto = 1.1 x 10~ 2 , 
with its position indicated by the vertical dotted line. 
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To understand our above observations, we 
gate probability distribution of the spin state, |c„ 

m\ip(t)}\ 2 . As shown in the insets of Fig. [21 we find 
that compared with the initial CSS, the maximally SSS 
at to has a very sharp probability distribution centered 
at the lowest spin projection, i.e., m = (for even N) 
or m = ±1/2 (for odd N). Such a sharp probability 
distribution of the SSS can be explained qualitatively by 
considering the familiar phase model [261 ] . By replacing 
J z -> P 4 = -idg, and J x -> (iVcos$)/2, with $ be- 
ing macroscopic phase difference between two condensate 
components, one obtain 



Ha, = -2k- 



d 2 

^$2 



SIrN 



■ cos<i>, 



(2) 



where we have taken fl(t) = to simulate quantum 
dynamics of the spin system before turning-off the field. 
The phase model Hamiltonian allows us to regard the 
spin system as a fictitious particle with effective mass 
(4k)" 1 subject to a pendulum potential. Moreover, the 
particle behaves as a pendulum rotating with frequency 
u> c s — v / 2kTTr]V in phase space ($, where mean 

phase difference $ = ($). As shown in Fig. UK a) and (b), 
starting from vertical distributed initial points ($ = tt 
and P$ = —j, —j + 1, one obtains the distribution 

elongated horizontally at the time t ~ T/4, where T = 
2ir/u> c ff is the period of the pendulum. It is reasonable 
to assume that the projection of the distribution along <!> 
will always be symmetrical to tt, e.g., <& = 2ir — $o and 
(Jx)o = jcos($ ), as Fig. [3(b). 
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FIG. 2: Schematic picture of the probability distribution in 
phase space ($, P<s>) for: (a) the initial CSS, (b) the SSS at 
t ~ T/4. The insets: the corresponding distribution |c m | 2 as 
a function of m (or P$) obtained numerically. 



Based upon the phase model, we explain intuitively 
why the obtained SSS has a sharp distribution with the 
lowest spin projection being occupied predominantly. In 
addition, we find that the sharp distribution accompanies 
with the maximal mean spin, i.e., (J x ): —j (J x )o a s 
$ : 7T — > $o, thus {d(J x )/dt)t Q = 0. On the other hand, 
from Heisenberg equation of J x we know the relation 



d(J x }/dt 



( J zJll 



JyJz 



-Atan(20 min ) [for A ^ 0, see 



the deduction of Eq. (6)] , so we further obtain 9 min = 
at to, as shown in Fig.[ljb). More important, we obtain 
analytical expression of the maximal-squeezing time 



T _ TT 
«*0 - K - - - y 2QrN , 



(3) 



which is valid for large N (> 10 3 ) and small coupling with 
k < il/f << Nk. By comparing with exact numerical 
results, we find that Eq. ([3|) gives accurate prediction of 
the maximal-squeezing time for flu near or larger than 
the optimal coupling. Law et al. have investigated the 
optimal coupling as a function of N based on a wide range 
of numerical simulations [l3j , from which we suppose that 
the optimal coupling obeys power rule £Ir/k ~ N 1 ^ 3 for 
the large N. Such as N = 10 3 , the optimal coupling is 
about Q R = 10.8k and Eq. © gives nt = 1.069 x 10" 2 , 
consistent with numerical result 1.104 x 10~ 2 . Note that 
Eq. ([3]) predicts the time scale for the sharp distribution 
of the SSS and the maximal (J x ) (also 9 m i n — 0), which, 
however, does not necessarily correspond to the maximal 
squeezing (see below). 

If the Josephson field is turned off at to, the spin sys- 
tem is governed only by the self-interaction Hamiltonian 
2k J 2 , so the distribution |c m (t)| is conserved while the 
relative phases among the spin projections are subject 
to change. When the SSS at the time to exhibits a 
very sharp distribution in to, as shown in the inset of 
Fig. [2Kb), the effect of relative phases induced by the 
self-interaction gives negligible influence to the squeez- 
ing. To show this, we suppose the SSS at to takes the 
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form 



Wo)) = 



e lv sin a 
for even N case, or 



(|i,l) + |j,-l)) + cosa|j,0), (4) 



m )) = — 3/2) - |i, -3/2» 
+£^(| il l/2)-|j I -l/2», 



(5) 



for the odd N, where tp and a represent the relative phase 
and the amplitude, respectively. Fig.[3]shows the squeez- 
ing parameter £ as a function of a and tp, where two 
distinct features are observed. Firstly, £ is minimized 
as a — ► 0], which implies that the maximal squeezing 
occurs as long as the SSS has a very sharp distribution 
with a large amplitude of the lowest spin projection or, 
equivalently, as the SSS approaches to the ground-state of 
2k J 2 . Secondly, £ looks insensitive to the relative phase 
ip for the SSS with a sharp distribution. Since the self- 
interaction only varies the relative phase ip with a fixed a, 
Fig. [3] explains qualitatively the storage of the maximal 
squeezing after to. 
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FIG. 3: The squeezing parameter £ as a function of a and ip 
(in units of 71") for (a) the even number case (N = 1000); (b) 
the odd number case (N = 1001), calculated by using Eqs. 
(|4]) and ([5]), respectively. 



The storage itself in our scheme does not depend on 
O/j, and the optimal coupling is chosen here to store the 
maximal squeezing £o- If is smaller than the optimal 
coupling, there exist two time scales: to for the vanishing 
# m i n , and tm for the maximal squeezing. Eq. ([3]) still 
works well to give the time scale of ^min = 0, but fails 
to predict that of the maximal squeezing. As shown in 
Fig. [5fa), for = 5k real maximal squeezing occurs at 
kto = 6.915 x 10~ 3 . From Fig. [4ja) we also find that 
the maximal (J x ) and (9 m j n = appears at the same time 
Kto = 1-687 x 10~ 2 , at which the probability distribution 
of the SSS is sharp enough [see the inset of Fig. HJb)]. 
As a result, a less squeezed variance with # m ; n = can 
be stored by turning-off the external field at the time to 
[see the red and the green lines of Fig. BJa)]. 
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FIG. 4: (Color online) Time evolution of £, 6 m i n and the mean 
spin (J x )/j for N = 10 3 , (a) Sl R = 5k, and (b) tl R = 50k. 
Black curves in (a) and (b) are the corresponding quantities 
for the constant-coupling case. Red (green) curves represents 
the squeezing £ (# m in) after turning off the field at time to- 
Inset: Probability distribution |c m | 2 as a function of m for 
the SSS at the times tq (the black line), and to with the blue 
line for (a) and the red line for (b). 



The validity of Eqs. ([4} and ([5]) can be tested by 
considering two exact solvable cases with N = 2 and 
N = 3. Taking the optimal coupling = k for 
the two-atom system and — 2k for the three- 

atom one, we obtain the maximally-squeezed states 
|V(t„)) = i{-l) n+1 e~ l ^\j = l,m = 0) i and 
W*»)> = ^e- 3i *W2 (| 3/2) 1/2) -13/2,-1/2)), re- 
spectively, where t n = (2n + l)ir/S for any integer n, 
and the level spacing S — 2y / Q 2 R + k 2 for N = 2 and 
S = 2^/fl' 2 R + 2k£1_r + 4k 2 for N = 3. Obviously, the spin 
state |-0(t„)) is the ground state of the self- interaction 
Hamiltonian 2k J 2 , which results in exactly constant £ 
with zero # m in by rapid turning-off the external field at 
the times t n . For large N, the SSS at the time to no longer 
corresponds to the ground state of 2k J 2 , but approaches 
to it compared with the initial CSS. Consequently, almost 
constant £ can be achieved after turning-off the external 
field. 



It is worth mentioning that the above qualitative ex- 
planation can not be applied to the SSS with a broad 
probability distribution. In fact, typical OAT-scheme [l| 
relies solely on the evolution of relative phases induced by 
the self-interaction, where the initial CSS shows a broad 
probability distribution. As shown by the red line of 
Fig. H{b), for a large coupling SIr — 50k, the squeezing 
parameter decreases slightly after turning-off the exter- 
nal coupling. This is because the SSS at to exhibits a 
relatively broad probability distribution [see the inset of 
Fig. 0Jb)]. From the green line of Fig.QJb), we find that 
the reduced variance is stored in the fixed direction with 
#min = 0. By comparing the numerical and the ana- 
lytical result of to for FIr — 50k, we also find that the 
numerical result 4.945 x 10~ 3 k -1 agrees very well with 
4.967 x 10" 3 k~ 1 estimated from Eq. ©. 

Before closing, we prove that the mean spin always 
appears in the x direction, and explain why we study 
the spin squeezing in the small-coupling regime. Note 
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that the linear combinations of the probability ampli- 
tudes pm\t) = c m (t) ± C- m (t) obey two closed sets of 
first-order ordinary differential equations. For even N, 
the fact that all p m (0) = results in pm\t) = 0, 
namely c_ m (t) = c rn (t). On the other hands, for odd 

N all pm\t) are zero, i.e., c_ m (i) = — c m (i). Since 
c-m(i) = ±c m (i), we obtain simple expressions: (J y ) = 
(J z ) = 0, and (J x ) 7^ 0, i.e., the mean spin is along 
the x direction. Consequently, the spin component nor- 
mal to the mean spin reads J n = J y sin 9 + J z cos 9. By 
minimizing the variance (AJ n ) 2 with respect to 9, we 



obtain the squeezing angle 9 

(A </n)min 



itan _1 (B/A) and 



= ±C - ±\M 2 + £ 2 , where A = (J 2 - J 2 ), 

£? = (JzJy + JyJ z ), and C = (J 2 + J 2 ). From Heisenberg 

equations of motion of the spin J a for a = x, y, z, one can 
obtain formal solutions for the constant coupling case: 

c = j{ 3 +i)-{Jl), a = -c+j(i-n R /K)-(n R /K){J x ), 

and B = -(2K)~ l d{J x )/dt. Note that for B = and 
A 7^ 0, the spin squeezing takes place along z axis (i.e. 
#min = 0) with the corresponding squeezing parameter 



H = i - (n R / K ) 



i + (X)o/j 



(6) 



where ( J^o is the maximum value of the mean spin. For 
an extremely strong coupling (Q, R >> kN), (J x )o —> —j 
and £o - * 1 so the squeezing becomes very weak. This 
is the reason why we discuss the spin squeezing in the 
small-coupling regime (k < fl R << Nk). 

Finally, we estimate several important parameters for 
experimental realization. Following Ref. [6j , we consider 
23 Na atoms in the hyperfine states \F — \,Mp = ±1) 
trapped in a spherically symmetric potential V a = V% = 
muj 2 r 2 /2. The self-interaction strength can be solved by 



applying the Thomas-Fermi approximation, yielding 
l5 2/5 hio a cS ( a ho N " ' 



14 aho \Na 



(7) 



where aho = y^h/mui being the harmonic oscillator 
length, and a c g = a aa + abb — 2a a b effective scatter- 
ing length. For 23 Na atoms, we take a aa = a^b = 
2a a b @ and Oeff = a>aa — 2.75nm [17], then the self- 
interaction strength n ~ 3.2448 x lO^^hw. For the 
case TV = 10 3 and Vl R = 10.8k, we have obtained 
to = 1.1041 x 10" 2 /(S" 1 k) = 34.03w-\ which corre- 
sponds to the maximal-squeezing time about 10.83 ms 
for lj = 2tt x 500Hz. 

In summary, we have investigated coherent control 
of spin squeezing by turning-off Josephson field at the 
maximal-squeezing time. We show that by applying the 
optimal coupling then turning off later, the maximal 
squeezing can be stored in the z axis for a long time, 
which can be explained in terms of the probability distri- 
bution of the SSS. For a sharp distribution with a large 
amplitude of the lowest spin projection, the effect of the 
self-interaction gives small contribution to the squeezing 
and its direction. We find the analytic expression of the 
maximal-squeezing time by considering the phase model. 
Our scheme for the storage of spin squeezing is quite ro- 
bust for a wide range of Josephson coupling. 

We thank Profs. OK. Kim, K. Nairn, C. P. Sun, 
W. M. Liu, S. Yi, and X. Wang for helpful discussions. 
This study was financially supported by Pusan National 
University in program Post-Doc 2006 and Korea Re- 
search Foundation Grant funded by the Korean Gov- 
ernment (MOEHRD, Basic Research Promotion Fund) 
(KRF-2006-312-C00543). 



[*] Present address: Department of physics, School of Sci- 
ence, Beijing Jiaotong University, Beijing 100044, China 

[t] Electronic address: swkim0412@pusan.ac.kr 

[1] M. Kitagawa and M. Ueda, Phys. Rev. A 47, 5138 (1993). 

[2] D. J. Wineland et al, Phys. Rev. A 46, R6797 (1992); 
ibid. 50, 67 (1994). 

[3] A. Kuzmich et al, Phys. Rev. Lett. 79, 4782 (1997). 

[4] J. Hald et al, Phys. Rev. Lett. 83, 1319 (1999). 

[5] J. M. Geremia et a/., Science 304, 270 (2004). 

[6] A. S0rensen et al, Nature (London) 409, 63 (2001). 

[7] K. Helmerson and L. You, Phys. Rev. Lett. 87, 170402 
(2001); M. Zhang et al, Phys. Rev. A 68, 043622 (2003). 

[8] X. Wang and B. C. Sanders, Phys. Rev. A 68, 012101 
(2003). 

[9] J. K. Korbicz et al, Phys. Rev. Lett. 95, 120502 (2005). 
[10] S. Yi and H. Pu, Phys. Rev. A 73, 023602 (2006). 
[11] U. V. Poulsen and K. Molmer, Phys. Rev. A 64, 013616 
(2001). 

[12] M. Takeuchi et al, Phys. Rev. Lett. 94, 023003 (2005). 
[13] C.K. Law et al, Phys. Rev. A 63, 055601 (2001). 



[14] S. Raghavan et al, Opt. Commu. 188, 149 (2001). 

[15] D. Jaksch et al, Phys. Rev. A 65, 033625 (2002). 

[16] D. S. Hall et al, Phys. Rev. Lett. 81, 1539 (1998); idib 

81, 1543 (1998). 
[17] J. Stenger et al, Nature 396, 345 (1998). 
[18] G. J. Milburn et al, Phys. Rev. A 55, 4318 (1997). 
[19] A. Smerzi et al, Phys. Rev. Lett. 79, 4950 (1997). 
[20] J.I. Cirac et al, Phys. Rev. A 57, 1208 (1998). 
[21] M.J. Steel and M.J. Collett, Phys. Rev. A 57, 2920 

(1998) . 

[22] E. M. Wright et al, Phys. Rev. Lett. 77, 2158 (1996). 
[23] D. Gordon and C. M. Savage, Phys. Rev. A 59, 4623 

(1999) . 

[24] A. Sinatra and Y. Castin, Eur. Phys. J. D 8, 319 (2000). 
[25] G. S. Agarwal and R. R. Puri, Phys. Rev. A 39, 2969 

(1989); G. R. Jin and W. M. Liu, ibid. 70, 013803 (2004). 
[26] D. Jaksch et al, Phys. Rev. Lett. 86, 4733 (2001); C. 

Menotti et al, Phys. Rev. A 63, 023601 (2001); A. 

Micheli et al, ibid. 67, 013607 (2003). 



